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Abstract 

Recently, a semiempirical method for alloys based on equivalent crystal theory was 
introduced (Phys. Rev. B45, 043(1992)). The method successfully predicts the con- 
centration dependence of the heat of formation and lattice parameter of binary alloys. 

In this report, a study of the parameters of the method is presented, along with new 
results for 7 Fe-Pd and 7 Fe-Ni alloys. 
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Recently, a new method for calculating alloy properties was introduced by Bozzolo, 
Ferrante and Smith (BFS) 1 . The method is based on the ideas of equivalent crystal theory 
(ECT) for defect formation energies in elemental solids 2 and uses only pure metal and two 
alloy properties as input data. In our previous work, we applied the method to the study 
of the heat of formation as well as the concentration dependence of the lattice parameter 
of several fee binary alloys 1,3 . More recent applications of BFS deal with the study of 
the energetics of bcc alloys 4 and theoretical modelling of the atomic force microscope for 
bimetallic 'u-. sample interactions 5 . Also, as a consequence of the ideas underlying the 
foundation of BFS, we were able to derive a new set of sum rules which allow for the direct 
calculation of bulk alloy properties from pure component properties 6 . 

The simplicity of the method relies on the basic assumption that the energy of each 
non-equivalent atom is described as a superposition of two separate contributions: a strain 
energy that deals with the structural changes and a chemical energy that takes into account 
the changes in chemical composition. The calculation of the first term, the strain energy, is 
a straightforward application of equivalent crystal theory for pure elements: it is computed 
as if the neighbors of a given atom were of the same atomic species. Thus, no informa- 
tion for the alloy is needed except for neighbor locations. The calculation of the chemical 
energy follows an ECT-like format with the introduction of two additional ‘perturbative’ 
parameters which mimick the interaction between atoms of different atomic species in the 
overlap region and account for the changes in the electron density due to the presence of the 
other atomic species 1 . For binary alloys A-B, two such parameters are needed, A ab and 
Ab. 4- which are obtained from two experimental (or theoretical) alloy properties as input 1 . 
Our previous work 1 3-6 only cited the application of BFS to several problems in order to 
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verify its ability to predict properties of interest. We now give a more detailed discussion 
concerning the properties of A ab and Aba and provide a set of values for general use. 

In our previous work 1,3-5 , Aab and Aba were determined by requiring that the pre- 
dicted heats of solution in the dilute limit reproduced the experimental values. These 
quantities were computed by making use of the cluster expansion method 7 which relates 
properties of ordered compounds to those of disordered. We used it to predict the heat of 
formation for disordered alloys. Although this scheme has been very successful in a large 
number of applications 8-10 , there is still some uncertainty regarding the truncation criteria 
used in the expansion 11 . This arbitrariness results in fluctuations of the predicted values 
of Aab and Aba ■ In some cases, the choice of a basis set on which the expansion is based 
could lead to theoretical predictions that differ greatly from the experimental ones, even if 
the physical quantity used as input (i.e., the heat of solution) is accurately reproduced. We 
will give an example of this along with providing a set of parameters Aab and Aba, used 
in our previous work 1,3-5 , and comment on the sensitivity of these parameters to the choice 
of truncation scheme used in the cluster expansion. 

Table 1 displays the pure element parameters needed for computing the strain and 
chemical energy contributions 1-2 and Table 2 shows the parameters A ab and Aba for 
several fee binary alloys of the elements listed in Table 1, obtained by using the tetrahedron 
approximation 8 (i.e., all the clusters contain only nearest-neighbors). In all cases we used 
the experimental heats of solution in the dilute limit as input, as this quantity is readily 
available for most alloys 12 . Previous applications of this nearest-neighbor model to total 
energy of bulk alloys can be found in the work of Connolly and Williams 8 , Terakura et al 9 
and Takizawa et al 10 . The ordered structures included in order to determine the many-body 
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potentials in the cluster expansion are the fee L\q and XI 2 structures. It was found that 
this basis set gives good results for all the alloys listed in Table 2. In tables 3 and 4 we 
display the corresponding parameters for bcc pure elements and binary alloys respectively. 

To the fee binary alloys discussed in Ref. 1 , we added two new cases: 7 Fe-Pd and 
7 Fe-Ni. The qualitative and quantitative agreement are good in both cases. Fig. 1 shows 
the heat of formation as a function of composition for these two disordered alloys. 

In table 5 we address the issue of stability of the parameters A with respect to the 
choice of basis set of ordered compounds in the cluster expansion for a Cr-V alloy 4 . The 
different choices are related to two possible ordered structures at 50 % composition (B 2 
and 532) and the corresponding pair multisite correlation functions ( £2 and { 3 ). We will 
denote the cases studied as follows: (i) 532-& includes the 532 ordered structure and the 
£2 correlation function (nearest-neighbor pair); (ii) 532-£j; 52-£> and (iv) B2 + 532-^2 + £3 
(see Ref. 4 for details). We note that the ambiguity in the values of such parameters is a 
consequence of the cluster expansion method and basis set chosen. Otherwise, we would 
expect to be able to define a unique set of parameters A for different alloys if no approximate 
schemes (for example, ab initio calculations) were involved in the theoretical determination 
of the property of choice. Figure 2 displays some results for the heat of formation of the 
Cr-V alloy (see Ref. 4 for a more detailed discussion). The several curves shown correspond 
to different choices of the basis set chosen, although the same experimental input is used in 
all cases. Obviously, the extreme difference in behavior of the heat of formation is purely 
a consequence of the basis set used and not of the formalism presented here. Actually, the 
parameters obtained for different basis sets are fairly stable: the parameters A play the role 
of perturbations on the EC'T parameter a (see table 1 and 3), giving a net change in a of 
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Element 

Cohesive 

Energy 

Lattice 

Constant 

P 

Ifl 

Q 

Hi 

A1 

3.34 

4.05 


0.336 

2.105 


Cti 

3.50 

3.615 


0.272 

2.935 

0.765 

Ag 

2.96 

4.086 


0.269 

3.337 

0.756 

Au 

3.78 

4.078 

10 

0.236 

4.339 

0.663 

Ni 

4.435 

3.524 

6 

0.270 

3.015 

0.759 

Ir 

6.94 

3.84 

10 

0.235 

4.408 

0.661 

Pd 

3.94 

3.89 

8 

0.237 

3.612 

0.666 

Pt 

5.85 

3.92 

10 

0.237 

4.535 

mEk 

Fe 

4.27 

3.57 

6 

0.279 

2.963 

ESI 


Table 1: Experimental input: Cohesive energy (in eV), lattice parameter (in A), ECT 
parameters: p. / (in A), a (in A -1 ) and A (in A) for several fee elements. 


the order of 0.01%. Thus in spite of the ability to obtain quantitatively accurate predictions 
in many cases it would be desirable to base the values of A on first-principles calculations. 
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Table 2: Parameters Aab and Aba for several fee binary alloys. 




Table 3: Experimental input: Cohesive energy (in eV). lattice parameter (in A). ECT 
parameters: p, / (in A), o (in A -1 ) and A (in A) for several bee elements. 
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A-B 

&AB 

Aba 

Cr-Fe 

0.0465 

0.0285 

Cr-Mo 

-0.02447 

-0.0090 

Cr-V 

-0.0246 

-0.0232 

Fe-V 

0.0998 

-0.07168 


Table 4: Parameters Aab and A ba for several bcc binary alloys. 


Basis 

Aab 

A BA | 


EMM Jj 

BiAlnil 



EE il 



-0.0222 


-0.0237 

-0.0225 


Table 5: Parameters Aab and Aba for different basis sets used in the cluster expansion 
(see Ref. 4) for Cr-V. 
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(a) Fe-Ni. 



(b) Fe-Pd. 

Figure 1 .—Experimental (solid squares) and BFS 
(solid line) values of the heat of formation (in 
eV/atom) versus concentration. 



Figure 2.— Experimental (soft) squares) and BFS 
(solid and dashed lines) values of the heat of 
formation On eV/atom) versus concentration for 
Cr-V. The different curves are obtained with 
different choices of the set of ordered structures 
considered. See Ref. 4 for a description of the 
different choices. 
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